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FOREWORD 


The  research  was  performed  under  Program  Element  62601 F,  Project  8809, 

Task  04. 

Inclusive  dates  of  research  were  1  October  1972  through  1  February  1974. 

The  report  was  submitted  6  May  1974  by  the  Air  Force  Weapons  Laboratory  Project 
Officer,  Major  Arthur  R.  Geldbach  (SAB). 

The  computer  subroutine  XAREA  has  been  developed  within  the  Air  Force 
Weapons  Laboratory  to  function  as  a  part  of  the  advanced  computer  model  QUANT0-- 
A  CODE  TO  OPTIMIZE  WEAPON  ALLOCATIONS  (AFWL-TR-73-242).  The  approach  used  in 
XAREA  for  providing  the  numerator  portion  of  the  P|<  based  on  overlapping  weapon- 
aircraft  area  has  been  reviewed  by  the  Systems  Analysis  Section  of  the  Battle 
Environments  Branch,  AFWL/SA,  and  is  considered  appropriate  for  use  in  activities 
relating  to  aircraft  flush  methods. 

The  basic  model  was  developed  upon  an  identification  of  need  by  Major 
Richard  Conway.  Assistance  in  the  formulation  of  the  unique  branching  mechanism 
was  provided  by  Captain  Karl  T.  Benson  who,  along  with  Major  Arthur  Geldbach, 
assured  the  completion  of  this  report  through  their  technical  editing  and 
experience. 

This  technical  report  has  been  reviewed  and  is  approved. 
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ABSTRACT 

(Distribution  Limitation  Statement  A) 


When  aircraft  are  flushing  from  a  base  under  attack,  the  portion  of  the  force 
that  will  escape  to  carry  out  the  mission  must  be  known  to  a  relatively  high 
level  of  confidence.  This  involves  the  determination  of  the  probability  of 
kill  associated  with  the  base  force.  One  major  component  of  this  term  is  the 
geometrical  intersection  of  the  weapon's  effective  kill  volume  with  the  region 
containing  all  aircraft  at  any  instant  in  time.  Subroutine  XAREA  calculates  a 
numerical  value  for  this  intersection  in  terms  of  area  according  to  the  follow¬ 
ing  assumptions:  (1)  all  interaction  takes  place  in  a  plane,  parallel  to  the 
ground  and  elevated  to  the  average  aircraft  altitude;  (2)  the  figure  within 
which  all  aircraft  are  located  is  an  annulus,  bounded  by  the  maximum  effective 
ground  distance  reached  by  the  first  aircraft  off  the  base  and  the  minimum 
distance  achieved  by  the  last  aircraft  in  the  flush  sequence;  (3)  the  effective 
kill  area  of  a  weapon  with  an  assumed  spherical  destruction  pattern  is  the 
circle  cut  by  tiie  plane  through  the  sphere.  Although  the  intended  use  of  this 
computer  function  is  highly  specialized,  it  is  applicable  under  any  conditions 
requiring  the  area  of  intersection  of  an  annulus  with  a  circle. 
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ABBREVIATIONS  AND  SYMBOLS 


Pt  Tiie  probability  that  an  aircraft  under  attack  will  be  destroyed  or 

K  "killed. " 

XAREA  The  intersection  area  of  a  specified  annulus  with  a  circle. 

AH  The  outer  radius  of  an  annulus. 

AL  The  inner  radius  of  an  annulus. 


XL 

XD 

XMIN 

XMAX 

ALANG 

AHANG 

XLANG 

OPANG 


The  radius  of  a  circle  which  intersects  an  annulus  to  produce  XAREA. 

The  displacement  distance  of  the  center  of  the  annulus  from  the 
center  of  the  intersecting  circle. 

In  a  rectangular  coordinate  system  oriented  about  the  center  of  an 

annulus  with  its  independent  axis  passing  through  the  center  of  the 

intersecting  circle,  this  is  the  horizontal  (x)  coordinate  of  the 
points  generated  from  the  intersection  of  the  circle  with  the  inner 
boundary  of  the  annulus. 

In  a  rectangular  coordinate  system  oriented  about  the  center  of  an 

annulus  with  its  independent  axis  passing  through  the  center  of  the 

intersecting  circle,  this  is  the  horizontal  (x)  coordinate  of  the 
points  generated  from  the  intersection  of  the  circle  with  the  outer 
boundary  of  the  annulus. 

The  angle  formed  between  the  line  segments  from  (1)  the  center  of 
the  annulus  to  the  center  of  the  intersecting  circle  and  (2)  the 
center  of  the  annulus  to  one  of  the  inner  annulus  intersection 
points. 

The  angle  formed  between  the  line  segments  from  (1)  the  center  of 
the  annulus  to  the  center  of  the  intersecting  circle  and  (2)  the 
center  of  the  annulus  to  one  of  the  inner  annulus  intersection 
points. 

The  angle  formed  between  the  line  segments  from  (1)  the  center  of 
the  intersecting  circle  to  the  center  of  the  annulus  and  (2)  the 
center  of  the  intersecting  circle  to  the  annulus  intersection  point 
(either  inner  or  outer  depending  on  where  it  is  applied). 

The  complementary  angle  opposite  one  of  the  above  angles  that  is 
greater  tren  90  degrees  when  such  a  condition  exists. 
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SECTION  I 
INTRODUCTION 

The  Air  Force  Weapons  Laboratory  (AFWL)  has  investigated  weapon  distribu¬ 
tions  and  their  relationship  to  kill  probability,  P^,  wherein  attention  was 
focused  on  the  effectiveness  of  a  weapon  with  a  spherical  destruction  pattern 
detonated  in  a  field  of  aircraft  flying  from  an  airbase  under  attack.  These 
aircraft  flushed  single  file  and  were  commanded  to  fly  at  different  headings 
by  executing  one  single  turn  up  to  180  degrees  in  any  direction  after  achieving 
minimum  required  turn  altitude  and  velocity.  In  this  manner,  they  were  con¬ 
strained  to  a  region  bounded  by  the  maximum  distance  reached  by  the  first  air¬ 
craft  off  the  base  and  the  minimum  distance  achieved  by  the  last  aircraft  in 
the  flush  sequence.  By  examining  all  possible  positions  the  first  aircraft 
could  reach  in  a  given  time  interval  for  any  turn  angle,  a  two-dimensional 
figure  can  be  generated  representing  the  ground  surface  a  .a  over  which  all 
aircraft  are  located.  As  flight  time  progresses,  this  figure  approximates  a 
circle  with  the  center,  or  centroid,  located  at  some  positive  distance  from 
the  point  of  brake  release  in  the  direction  of  takeoff.  Given  enough  time,  the 
final  aircraft  in  the  flush  sequence  will  escape  from  t'  base,  leaving  a  void 
within  which  no  aircraft  are  located.  As  flight  time  continues,  this  void 
region  approximates  a  circle  concentric  with  the  maximum  distance  (i.e.,  first 
aircraft's)  circle,  but  of  smaller  radius.  Thus,  after  a  period  of  elapsed 
time,  the  region  containing  all  aircraft,  referred  to  as  the  annulus,  is  the 
area  inside  a  large  circle  and  the  area  outside  a  smaller  circle,  both  circles 
being  centered  about  the  same  point.  For  simplicity,  this  annular  region  is 
assumed  to  be  contained  in  a  horizontal  plane  at  some  aircraft  altitude.  If  a 
weapon  with  an  effective  kill  volume  is  detonated  at  random  in  space,  there 
will  be  no  chance  of  destroying  any  of  the  aircraft  unless  the  weapon  volume 
intersects  the  region  where  the  planes  art  located.  Should  this  happen,  then 
the  kill  probability  can  be  computed  basea  on  the  overlapping  area.  Assuming 
the  weapon's  lethal  volume  is  spherical,  then  the  lethal  region  in  the  horizon¬ 
tal  plane  containing  the  aircraft  is  a  circle  (figure  1).  If  the  aircraft  are 
uniformly  distributed  w J L!n i>  the  annular  region,  the  probability  of  killing  an 
aircraft  is  the  quotient  of  (1)  the  area  common  to  both  the  weapon's  circular 
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lethal  area  and  the  aircrafts'  annular  region,  divided  by  (2)  the  area  of  the 
annulus  containing  all  aircraft.  Subroutine  XAREA  has  been  created  to  compute 
the  numerator  of  this  quotient,  for  all  combinations  of  values  of  the  two  radii 
of  the  annulus,  the  radius  of  the  lethal  circle  of  the  weapon,  and  che  offset 
distance  of  the  weapon  circle  center  form  the  center  of  the  annulus. 


Figure  1.  Horizontal  Annulus  Plane  Intersecting  Lethal  Weapon  Volume 
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SECTION  II 

MATHEMATICAL  DEVELOPMENT 


Subroutine  XAREA  computes  the  area  common  to  a  circle  and  an  annulus.  The 
value  returned  to  the  main  program  is  the  area  of  intersection  of  the  annular 
region  of  outer  radius  AH  and  inner  radius  AL,  with  a  circular  weapon  d^ea  of 
lethal  radius  XL  displaced  a  distance  XD  from  the  centroid.  The  centers  of 
the  circles  defining  the  annulus  are  concentric  with  the  centroid  of  the  air¬ 
craft  area.  It  should  be  noted  that  in  all  cases  AH  >  AL  and  all  measurements 
are  taken  to  be  nonnegative  with  respect  to  the  centroid,  in  the  same  unit 
denominations. 

The  basic  formulas  used  for  the  computation  of  the  area  of  intersection 
involve  repeated  application  of  the  rules  governing  the  area  of  circles  and 
segments  of  circles.  The  equations  for  these  areas  are 


A  (circle)  =  ttR2 

(1) 

A  (segment)  =  1/2R2(0  -  sin  6) 

(2) 

where  R  is  the  radius  of  the  circle  and  9  is  the  central  angle  (in  radians)  of 
the  segment.  The  circles  are  taken  two  at  a  time,  and  a  triangle  is  constructed 
with  vertices  at  the  center  of  each  circle  and  one  of  the  two  points  where  the 
circles  intersect  each  other  as  shown  in  figure  2.  In  this  example  the  three 
sides  of  the  triangle  are  indicated  as  XD  (distance  between  centers),  XL  (dis¬ 
tance  from  weapon  detonation  point  to  circular  intersection  point),  and  AL 
(distance  from  centroid  to  circular  intersection  point).  Any  angle  in  a 
triangle  may  be  found  from  the  formula 


Angle  =  arc  cos 


■ 

(adjacent  side)2  +  (other  adjacent  side)2  -  (oppusite  side)2 
2  •  (adjacent  side)  •  (other  adjacent  side) 


(3) 


By  determining  all  the  angles  of  the  triangle  AL-XL-XD  and  finding  the  area  of 
all  segments  involved,  a  quick  examination  of  the  orientation  of  the  intersection 
will  reveal  which  segments  should  be  added  and  which  should  be  subtracted  to 
obtain  the  correct  solution. 
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Assume  that  by  using  the  above  formulas  the  areas  have  been  Pound  for 
segments  A  and  8  in  the  orientation  shown  in  figure  2.  The  area  of  intersec¬ 
tion  of  the  aircraft  ar.nulus  with  the  weapon's  lethal  circle  will  clearly  be 
the  area  of  segment  A  minus  the  area  of  segment  B.  These  areas  arise  from  the 
application  of  formula  (2),  where  XL  represents  R  when  calculating  segment  A, 
and  AL  represents  R  for  the  determination  of  segment  B.  Theta  is  angle  B  or  a 
in  figure  2,  respectively,  and  is  found  from  a  direct  application  of  formula 
(3)  using  the  three  triangle  sides  AL,  XL,  and  XD.  In  the  program  itself  these 
angles  are  named  for  quick  identification.  ALANG  (AL  angle)  indicates  the 
angle  formed  with  AL  and  XD  as  adjacent  sides  ( a  in  this  case).  XLANG  (XL 
angle)  is  similar  but  has  XL  and  XD  as  adjacent  sides  (B  in  figure  2).  Both 
of  these  angles  have  XD  as  a  common  side.  Note  that  the  area  of  segment  A 
does  not  follow  directly  from  formula  (2).  The  segment  area  found  from  using 
angle  8  and  triangle  side  XL  must  be  subtracted  from  the  area  of  the  lethal 
circle  to  obtain  the  actual  area  of  segment  A.  This  subtraction  is  not 
necessary  for  determining  the  area  of  segment  B. 

In  this  first  example  the  angles  of  triangle  AL-XL-XD  are  all  acute 
angles.  It  is  possible,  however,  for  one  of  these  angles  to  be  larger  than 
90  degrees  under  certain  geometrical  conditions.  The  second  example  depicting 
this  situation  is  identical  to  the  first  one  with  one  alteration.  The  maximum 
distance  the  first  aircraft  off  the  base  could  achieve  (AH)  has  been  decreased 
so  that  the  lethal  circle  entirely  avers  a  portion  of  the  annulus  (figure  3). 
Although  the  computation  of  the  area  of  segments  A  and  B  is  identical  to  the 
previous  example,  two  new  segments,  l  and  D,  must  be  analyzed  to  determine  the 
true  intersectional  area,  for  a  portion  of  segment  A  furthest  from  the  centroid 
extends  beyond  the  intersection.  Again,  formula  (2)  is  used  to  find  these  new 
segment  areas.  This  time  the  radius  R  is  represented  by  AH  for  segment  C  and 
XL  for  segment  D.  The  angles  m  and  represent  6  in  this  equation,  where  ip  is 
derived  by  calculating  the  complimentary  angle,  4>,  from  formula  (3)  and  sub¬ 
tracting  the  result  from  180  degrees.  In  the  program,  angle  \i  is  denoted  as 
AHANG  (AH  angle),  ip  is  again  XLANG,  and  <p  is  called  OPANG  (opposite  angle) 
since  it  is  complementary  to  angle  ip. 
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SECTION  III 
THE  PROGRAM 

There  are  several  conditions  surrounding  the  intersection  of  a  circle  with 
an  annulus.  It  is  possible  for  the  lethal  circle  to  be  smaller  than  the  inner 
annulus  circle,  larger  than  the  outer  annulus  circle,  or  between  ty  *  two  in 
relative  size.  Within  this  context,  the  entire  lethal  circle  could  be  small 
enough  to  fit  within  the  annulus,  or  large  enough  to  overlap  both  the  inner 
and  outer  dimensions  of  the  annulus  if  properly  oriented  between  the  two.  The 
same  basic  formulas  are  applicable  in  each  case,  but  different  approaches  are 
required  in  summing  the  individual  segment  areas  together  properly  to  determine 
the  actual  area  of  intersection.  It  must  be  known,  for  example,  if  the  lethal 
circle  actually  intersects  the  annulus.  If  it  does,  it  should  be  determined  if 
the  intersection  is  with  the  inner  annulus  circle,  the  outer  annulus  circle, 
or  both. 

In  all,  there  are  seven  different  categories  of  tests  required  to  isolate 
the  geometrical  orientation  of  a  particular  situation  and  define  the  proper 
mathematical  scheme  required  for  intersection  area  calculation.  Within  the 
program  five  of  these  categories  have  only  one  possible  solution  to  a  given 
series  of  questions  presented  as  "IF"  statements.  A  numerical  value  is 
assigned  to  each  of  the  questions,  ranging  from  one  up  to  the  total  number  of 
questions  in  the  category.  Only  the  numerical  value  for  the  correct  solution 
to  one  question  in  each  category  is  carried  through  the  program,  adjusted  by 
a  factor  of  10  to  distinguish  the  categories  from  each  other.  At  the  completion 
of  the  questioning  period  the  five  resultant  numbers  are  summed  up.  Each  digit 
of  this  final  solution  represents  the  correct  answer  from  each  category.  This 
number  is  then  compared  against  an  array  of  possible  number  combinations  until 
the  mating  value  is  located.  At  that  point  the  control  transfers  to  one  of 
the  solution  blocks  for  actual  intersection  area  computation  through  the  use 
of  a  computed  GO  TO  statement.  If  no  matches  have  occurred,  the  subroutine 
will  return  an  area  of  -1.0  as  an  error  indicator  for  invalid  input  data  since 
all  proper  solutions  are  nonnegative. 
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Appendix  I  has  been  provided  to  help  trace  through  XAREA  and  find  the 
section  of  code  involved  in  calculating  the  intersection  area  for  any  specific 
example.  By  answering  the  questions  in  each  test  section  using  numerical 
values  for  the  four  input  variables,  the  corresponding  statement  label  of  the 
start  of  the  pertinent  solution  coding  in  XAREA  can  be  identified.  By  following 
through  this  coding  block  in  the  program  listing  {appendix  III),  the  area  of 
intersection  of  the  annulus  with  the  lethal  circle  may  be  calculated  for  a 
specific  situation.  To  implement  subroutine  XAREA  on  a  computer  system,  a 
sample  driver  program  has  been  provided  in  appendix  II,  with  the  expected 
output  shown  in  appendix  IV. 
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SECTION  IV 
SUMMARY 

Although  the  formulas  used  for  computing  the  area  of  circles  and  segments 
of  circles  are  simple,  there  are  several  geometrical  conditions  surrounding 
the  intersection  of  a  circle  with  an  annulus  that  require  much  more  detailed 
analysis.  In  light  of  these  complications,  and  with  a  strong  desire  to  be  able 
to  calculate  this  area  of  intersection  under  any  possible  condition,  a  separate 
computer  subroutine  was  developed  as  a  subtask  to  a  larger  Air  Force  study. 

This  function,  called  XAREA,  has  been  written  to  satisfy  a  definite  need 
recognized  by  the  Air  Froce  Weapons  Laboratory.  It  is,  however,  applicable 
to  any  situation  requiring  the  area  of  intersection  of  an  annulus  with  a  circle. 
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CASE  TEST  AND  EVALUATION 

The  following  tests  are  designed  to  direct  a  user  to  the  pertinent  portion 
of  the  subroutine  coding  where  the  actual  area  calculation  is  performed  for 
specific  case  situations.  The  input  parameters  required  before  beginning  the 
tests  are 

XL  The  lethal  circle  radius 

XD  The  offset  distance  of  the  lethal  circle  from  the  center  of  the 
annulus 

AL  The  radius  of  the  inner  annulus  circle 

AH  The  radius  of  the  outer  annulus  circle 

The  resultant  coding  location  may  be  found  in  appendix  III.  The  functional 
elements  ANGLE,  SEGMT1,  and  SEGMT2  are  defined  at  the  beginning  of  the  code. 

TEST  1 

Size  relationship  of  lethal  circle  to  annulus.  One  of  the  following 
conditions  will  hold: 

If  2XL  <  AH-AL  and  XL  <  AL,  write  down  10,000. 

If  2XL  <  AH-AL  and  AL  <  XL  <  AH,  write  down  20,000. 

If  2XL  >  AH-AL  and  XL  <  AL,  write  down  30,000. 

If  2XL  >  AH-AL  and  AL  <  XL  <  AH,  write  down  40,000. 

If  2XL  >  AH-AL  and  XL  >  AH,  write  down  50,000. 

TEST  2 

Relationship  of  left-most  edge  of  the  lethal  circle  to  the  annulus  for 
positive  values  of  XD  measured  to  the  right  of  the  center  of  the  annulus. 

One  of  the  following  conditions  will  hold: 

If  XD-XL  <  -AH,  add  1000  to  the  number  from  test  1. 

If  -AH  <  XD-XL  <  -AL,  add  2000  to  the  number  from  test  1. 

If  -AL  <  XD-XL  <  AL,  add  3000  to  the  number  from  test  1. 
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If  AL  <  XD-XL  <  AH,  add  4000  to  the  number  from  test  1. 

If  XD-XL  >  AH,  add  5000  to  the  number  from  test  1. 

TEST  3 

Relationship  of  right-most  edge  of  the  lethal  circle  io  the  annulus  for 
positive  values  of  XD  measured  to  the  right  of  the  center  of  the  annulus.  One 
of  the  following  conditions  will  hold: 

If  XO+XL  <  AL,  add  100  to  the  number  from  test  2. 

If  AL  <  XD+XL  <  AH,  add  200  to  the  number  from  test  2. 

If  XD+XL  >  AH,  add  300  to  the  number  from  test  2. 

TEST  4 

Determination  of  the  intersection  of  the  lethal  circle  with  the  outer 
annulus  circle.  An  intersection  exists  if  the  number  derived  from  tests  1,  2, 
and  3  matches  any  of  the  following: 

14300  24300  33300  42300 

34300  43300 

44300 

If  the  numbers  match,  then  calculate  XMAX,  an  indicator  for  obtuse  angles  in 
the  triangle  AH-XD-AL,  from 


XMAX  =  AH.2-  -  XD-2-  X.L.; 

2  •  XD 

If  the  numbers  do  not  match,  set  XMAX  equal  to  a  very  large  positive 
number  approximating  positive  infinity.  Due  to  computer  limitation 
+1,000,000,000  was  used.  Note  that  the  number  combinations  52300,  53300,  and 
54300  also  indicate  an  intersection  of  the  lethal  circle  with  the  outer  annulus 
circle;  however,  when  XL  >  AH,  this  calculation  of  XMAX  is  not  required. 

TEST  5 

Determination  of  the  intersection  of  the  lethal  circle  with  the  inner 
annulus  circle.  An  intersection  exists  if  the  number  derived  from  tests  1,2, 
and  3  matches  any  of  the  following: 
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13200  23200  33200  43200  53300 

33300  43300 

If  the  numbers  matcn,  then  calculate  XMIN,  an  indicator  for  obtuse  angles  in 
the  triangle  AL-XD-XL,  from 


XMIN  = 


AL2  +  XD2  -  XL2 
2  •  XD 


If  the  numbers  do  not  match,  set  XMIN  equal  to  a  very  large  negative 
number  approximating  negative  infinity.  Due  to  computer  limitation 
-1,000,000,000  was  used. 

TEST  6 

Branching  for  special  obtuse  angle  considerations  through  the  relationship 
of  the  center  of  the  lethal  circle  to  the  intersection  points  found  in  tests  4 
and  5.  One  of  the  following  conditions  will  hold: 

If  XD  <  XMIN,  add  '.u  to  the  number  from  test  3. 

If  XMIN  <  XD  <  XMAX,  add  20  to  the  number  from  test  3. 

If  XD  >  XMAX,  add  30  to  the  number  from  test  3. 

TEST  7 

Orientation  of  the  intersection  of  che  lethal  circle  with  the  inner  annulus 
circle  in  relation  to  the  center  of  the  annulus.  One  of  the  following  condi¬ 
tions  will  hold: 

If  XMIN  <  0,  add  1  to  the  previous  number. 

If  XMiN  >  0,  add  2  to  the  previous  number. 

TEST  EVALUATION 

The  five  digit  number  calculated  from  the  test  section  indicates  which  set 
of  mathematical  functions  coded  in  XAREA  are  used  to  calculate  the  actual 
intersection  area  for  the  geometrical  orientation  under  evaluation.  Match  this 
number  with  the  exact  same  value  listed  below  to  determine  the  label  location 
used  in  the  code  for  area  calculation.  If  none  of  the  numbers  match,  an  error 
has  been  made  in  either  the  input  data  or  the  test  section. 


12 


AFWL-TR-74-67 


Number 

Label 

Number 

Label 

Number 

Label 

13121 

3 

33121 

3 

43321 

16 

13212 

4 

33212 

4 

43322 

12 

13222 

5 

33222 

5 

43331 

17 

14221 

6 

33312 

11 

43332 

13 

14321 

7 

33322 

12 

44321 

7 

14331 

8 

33332 

13 

44331 

3 

15321 

3 

34321 

7 

45321 

3 

22221 

9 

34331 

8 

51321 

20 

23221 

10 

35321 

3 

52321 

15 

23222 

5 

42221 

9 

53312 

22 

24221 

6 

42321 

14 

53321 

21 

24321 

7 

42331 

15 

53322 

22 

24331 

8 

43221 

10 

54321 

3 

25321 

3 

43222 

5 

55321 

3 
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APPENDIX  II 

DRIVER  PROGRAM  LISTING 


PROGRAM  iJPlV'FK(lNPur»QOTPi]T»TAPEb=TNiPUT»  f  APE6=0UTPUT ) 

M  =  S 
N  =  b 

aL  lNFb=5'). 

1  kfc  AD  (M » )  aL»«L»AH,aD 
K  I-OWMAT  (4E  lij.bl 
IE (XL.LT.O.O)STOP 
*PlfF(N»3)  XL  ♦  AL  »  Ap 

3  f  OPM AT (iHltlOA* 3HXL  =  *F 10 .S.5X . 3HAL  =  ♦ F 1 O.S»5X «  IHAH= ♦ F 1 0 . S// 
I  1 6X  *  <?HXU«  1  OX  ♦  1  7HIN  fFPbcCT  IMG  A«F A ♦  1  ♦  2HPK/ ) 

H  APEA=xAPFA (XL»XU*Al» Ah) 

TAPtA  =  J.141S0*!654*<  AH**k!-AL**?> 

PK=Akt A/ T  APF A 
*PlTF(N»b)  XD»  APcA.PK 
d  F  OPM  AT  ( 10X  ,F10.5*  10X*H2.S*  10X*F10.7) 

XO=XQMAl*XL)/XLlNf-S 
IF  (  Xi). lT  .Ah»XL)  GO  TO  4 
bO  TO  1 
END 
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APPENDIX  III 
SUBFUNCTION  LISTING 


no  ooooooooooooooooononnooorjononoono 


FUNCTION  XAREA  ( XL ♦ AO » AL . AHI 
DIMENSION  NUMBER! 1 7 >  «  NUMBER (6) 

UAT A  P 1/3 « 1415926536/ 

UATA  (NUMBER! J) « J= 1 « 1 7 ) /2121 .222 1 . 232 1 . 2331 . 3121 , 3212 ♦ 3221 , 3222* 
1 1 2. 3321. 3322. 3331. 3332. *22 1.4321 .4331.5321/ 

UATA  (NUMBER <J).J= 1.61/51. 61. 7 I. 72. 8 1.91/ 


SUBROUTINE  X AREA  RETURNS  THE  INTERSECTION 
AREA  OF  AN  ANNULUS  WITH  A  CIRCLE. 


INPUT  VARIABLE  deschiption- 

XL  —  THE  RAOIUS  OF  THE  INTERSECTING  CIRCLE 
XD  —  THE  POSITIVE  OISTANCE  BETWEEN  THE 
CENTER  OF  THE  INTERSECTING  CIRCLE 

and  the  center  or  the  annulus 

AL  —  THE  OUTER  RADIUS  OF  THE  ANNULUS 
AH  —  THE  INNER  RftOIuS  OF  THE  ANNULUS 


WESTR ICT IONS- 

1.  ALL  NUMBERS  MUST  BE  REAL  AND  NON-NEGATIVE. 

2.  AH  MUST  faE  GREATER  THEN  AL. 

3.  ALL  NUMBERS  MUST  BE  IN  THE  SAME  UNITS. 


ROUTINE  DEVELOPED  BY  BILL  PEAV.  AFWL/SAB*  KAFB,  NEW  MEXICO 

12  OCTOBER  1972 
REVISED  24  JANUARY  1974  (WFP) 


DEFINED  FUNCTIONS— 

ANGLE (AA.BB. CC) =ACOS ( (AA**2»BB**2-CC**2 1/(2. 0*AA*6B) ) 
SEGMT1 (R.THETA)=R**2* (PI -(THETA-SIN (THETA)  1/2.0) 
SE6MT2(R»THETA)*(R'  *2*(THETA-SIN(THETA) 11/2.0 


SET  XMAX  MIN  REFERENCE  LEVELS 


XMAX=* 1 «E*9 
XMIN=-I.E*9 


WFP00I00 
wFPOOll!) 
WFP00120 
33WFP00130 
WFP00I40 
WFP00150 
WFP00160 
WFP001 70 
wFPOO 1 80 
WFP00I90 
WFP00200 
WFP00210 
WFP00220 
WFP00230 
WFP00240 
WFP00250 
WFP002S0 
WFP00270 
WFP00280 
AFP00290 
WFP00300 
JFP003I0 
WFP00320 
WFP00330 
WFP00340 
WFP00350 
WFP00360 
WFD00370 
WFf'00380 
WFP00390 
WFP00400 
WFP00410 
WFP00420 
WFP00430 
WFP00440 
WFP00450 
WFP00460 
WFP00470 
WFP00480 
WFP00490 
WFP00500 
WFP0051 0 
WFP00520 
WFP00530 
WFP00540 
WFP0055  J 
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c 

c 


c 


4 


5 


preliminary  branching 

IF  (AL.GF.AH)  on  I,)  2 
IF  (XL.GT.AH)  GO  (■'  i  8 

MAIN  Itsr  BLOC* 

test  fok  left  intersection 

A=XO-XL 

IF  (A.GF. (-AHJ .ANO.A.lF. < -AL ))  t4=2000 
IF  (A.GT  .  <-AL>  .ANO.A.lE.  (AL  )  )  14=3000 
IF  { A.GT  .AL.ANu.A.Ll  .AHi  14-4iKiO 
IF  (A.GT. AH)  U-SOuO 

TEST  FOR  R 1.  HT  INTERSECTION 

B=XD*XL 

IF  (ri.Li  .  Al  »  I  3=1 00 

IF  (H.GI . AL . ANG.ri.Lt.AH)  13=200 
IF  (H.GT.AM)  13=300 

TESTS  INVOLVING  XMAX  AND  XMIN 


IF 

IF 

IF 

IF 

IF 

IF 

IF 


(  I4.E).  10'tu.ANi).  1  3.GE  .200) 

(  I  3,1  ■  4  .  }l)0.  ANO.  14.LE  .4000! 

(  X  D  •  t  F  •  X  M  I N )  12=  1U 

(xo.ct . xm i n. and. xo.Lc xmaxi 

(XO.GT.XMAX)  1 2  = 30 
(XMIN. LF  .0.13)  11  =  1 

(XM1F  .01 .0.0)  11=2 


XMIN=(XD*°2* Al  »*2-XL**2>/ (2.0*XD) 
Xm4X=(X0**2*Ah»»('>-hL»*2)/(2.0*XD) 

12=20 


NLIM=  I4*I3*I?*II 


OETFRM1NE  COMPONENT  SUM 


CALCULATE  ROUTING  ARGUMENT 

I TST=1 
DO  1  J= 1 « 1 7 

IF  <NUM.GE.NUMHEH< j) )  ITST=ITST»I 
CONTI  NUT 


BRANCH  TO  appropriate  CODE  BLOCK 


GO  TO 


(2.3.9.14.1S.3.4.10.S.11.16. 12.17. 13.F.7.8.3). 


ITST 


CODE  BLOCK  SECT  1 0N-- 

XARF A  =- 1 .0 
RETURN 
XAREA=0 .0 
RETURN 

ALANg=ANGLF  (AL.XD.XL)»2.0 
OP ANG= ANGLE ( AL  *  XL ♦ XO) 

XLANG= (ALANG/2.0*OPANG) *2.0 
ASEGMT=‘>EGMT2(XL.XLAN(j) 

BSEGMT  =  SEGMT2 ( AL  t  Al  ANG  F 
XAREA=A  .EGMT-BStGMT 
RETURN 

AL ANG= ANGLE <AL*XO* XL) *2.0 
XLANG=ANGLE(XL.XO.AL)*2.0 
ASFGMT=SEGMT1 (XL.XLANG) 

BSEGMT=SEGMT2<AL.ALANG) 

X  ARt  A=A  "iE  GMT-BSE  GMT 
RETURN 


WFPOObBO 
WFR00S70 
WFP00S80 
WFP00B90 
WFP00600 
KFP00610 
WFP00620 
WFPCG630 
WFPO0640 
MFP006S0 
HFP00660 
WFP00670 
WFP006B0 
*FP00fc90 
*)FP00700 
WFPOOi  10 
WFPOO  i’20 
WFP00730 
WFP00740 
WFP00750 
WFPOO  760 
BFP00770 
^00780 
WFP00790 
WFP00800 
BFP00810 
WFP00820 
WFP00830 
WFP00840 
WF  P008S0 
HFP00860 
HFP00870 
WFP0&880 
IVFP00890 
WFP00900 
WFP00910 
*FP00920 
WFP00930 
WFP00940 
WFP00950 
WFP00960 
WFP00970 
WFP00980 
WFP00990 
WFP01000 
WFP01010 
WFP01020 
WFP01030 
WFP01040 
WFP01050 
WFP01060 
PO 1070 
*FP0 1080 
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6  XAREA=PI*XL**2 
RETURN 

7  AHANG=ANGLE(AH,XD*XL)*2.Q 
OP ANG= ANGLE ( AH,XL  *XD) 
XLANG=(AHANG/2.0*OPANG)*2.0 
ASEGMT=SEGMT1  (XL*XLAN(,) 
BSEGMT=S£GMT2<AH*AHANG) 

XAREA=ASEGMT ♦USEGM  f 
RETURN 

a  AHANG=ANGLE(AH.XO*XL)*2.0 

XL  ANG= ANGLE  (XL*  XU.  AN)  «2.0 
ASEGMT=SEGMT2(XL*XLANG) 
«SEGHT=SEGHT2(AH,/hANG) 

XAREA  =  ASEGMT  *bSEGMT 
RETURN 

9  XAREA-PI*(XL**2-AL**2) 

RETURN 

10  XLANG=ANGLE(XL.XO> AL)*2.0 
OP ANG= ANGLE (XL.AL*XD) 
ALANG=(XLANG/2.0*OPANG)*2.0 
ASEGMT=SEGMT2<AL»AlANG) 
BSEGHT=SEGMT2(XL.XlANG) 
XAREA=PI*(XL**2-AL**2) ♦ASEGMT-BSFGMT 
RETURN 

11  ALANG=ANGLE(AL»XO*XL)*2.0 
OPANG=ANGLE ( AL . XL . XU) 
XLANG*(ALANG/2.0*OPANG)*2.0 
ASEGMT=SEGMT2(XL.XLANG) 

8SEGMT=Sf GMT2(AL.ALANG) 

AH ANG* ANGLE (AHtXO*XL)*2«Q 

OPANG=ANGLE ( AH, XL , XO) 

XLANG=(AHANG/2.0*OPANG>*2.0 

CSEGMT*SEGMT2(XL,aLANG) 

0SEGMT=SEGMT2(AH,AhANG) 

XAREA=ASEGMT*USEGMT-BSEGMT-CSFGMT 

RETURN 

12  ALANG=ANGLE  (AL*XO,XU*2.0 
aLANG=ANGLE ( XL *X0» AL ) *2. 0 
ASEGMT=SEGMT1 (XL*XLANG) 
8SEGMT=SEGMT2(AL»AlANG) 

AHANG=ANGLE ( AH»XD*XL>*2.0 

opang=angle ( AH , XL  » XO ) 

XLANG= (AHANG/2»0*OPANG) *2.0 

CSEGHT=SFGMT2(XL»XlANG) 

0SEGMT=S£GMT2(AH,AHANG> 

xarea=asegmt*osegmt-8segmt-csegmt 

RETURN 


WFP01090 
WFPO 1100 
WFP01110 
WFP01120 
WFP01 130 
WFP01140 
WFP01150 
WFP01160 
WFP01170 
WFP01180 
WFP01190 
WFP01200 
WFP01210 
WFP01220 
WFP0123Q 
WFP01240 
WFP01250 
WFP01260 
WFP01270 
WFP01280 
WFP01290 
WFP01300 
WFP01310 
WFP01320 
WFP01330 
WFP01340 
WFP01350 
WFP01360 
WFP01370 
WFP01380 
WFP01390 
WFP01400 
WFP01410 
WFPO 1420 
WFP01430 
MFP01440 
WFP01450 
WFP01460 
WFP01470 
•TPO 1480 
WFP01490 
WFPO 1500 
WFP01510 
WFP01520 
WFP0153C 
WFPO 1540 
WFP01550 
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13  ALANG=ANvLE(AL.XD.XL>*2.0 
Xl_ANG=ANGl  E  (XL.XU.AL>*2.0 
ASEGMT=SEGMT?  (XL.XLANG) 

BSEGM I =Sf  0MT2 ( AL  * AL  ANG ) 

AWANG=ANGlE (AH,XD.XL)*2.0 
XL ANG= ANGLE (XL.XO.AH) *2.0 
CSEGMT=Sf  i,MT2(XL.XlANG> 

DSEGMT  =SE  GMT2 (  AH*  Am  ANG ) 
XAREArCSEGMT*DSEGMT~ASEGMT~8SFGMT 
RE  TURN- 

14  amANG=ANGLE (AH,XD.xL>*2.0 
OR ANG = ANGLE (Art.XL.XO) 

XLANG=(AHAnG/2.0*OPANG)*2.0 
ASEGMT  =SEGHT 1 (XL.XLANG) 

HSEGMT=SEGMT  2 ( AH*  AH ANG 1 

XAREA=ASEGHT»BSEGMT-PIftAL**2 

RETURN 

lb  AHANG=ANGl  t(AH,Xl),xL>*2.0 

XLANG  =  ANGl.E  ( XL  «  XU.  AH)  *2 . 0 
ASEGMT  =  SEi»MT?(  XL  *XL  ANG) 

HSF.GMT  =  SEGMT2(AH,AhANG> 

X ARE A= ASEGMT ♦bSEGMT“P !*AL,#2 
RETURN 

16  XLANG=ANGL£  (  XL  »  XI).  AL )  *2  •  0 
OPANG=ANC-LE(XL.AL.xO> 

ALANG=(XLANG/2.0*OPANG)*2.0 

ASEGMT=SEGMT?(ALiALANG> 

BSE GMT=SEGMT2( XL.XLANG) 

AH ANG= ANGLE ( AH,  XI).  XL )  *2.0 
OPANG=ANGLE (AM. XL.XO) 

XLANG=( AhANG/2.0*OPANG)*2.0 
CSEGMT  =SF GMT2 (XL.XLANG) 

DSEGMT=SEGMT2(AH,AHANG> 

XAREA=Pl® (XL**2-AL**2> ♦ASEGMT.DSEGMT-bSEGMT-CSEGMT 
RETURN 

17  XL  ANG= ANGLE  (  XL  ,  XI) »  AL)  *2 .0 
OPANG=ANGLE(XL.AL,XD> 

ALANG= (XLANG/2.0*GPANG)*2.0 
ASEGMT=SEGMT2( al.alang) 

8SEGMT=SEGMT2 (XL.XLANG) 

AH ANG= ANGIE (Ah, XU. XL) *2.0 
XL  ANG=  ANGLE  (XL.  XI),  AH )  *2.0 
CSEGMT=SEGMT.?  (XL  )  XL ANG) 

OSEGMT  =  SEGMTi>  (AH.AHANG) 

XAREA=ASF GMT .CSEGMT «OSEGMT-BSEGMT-PI*AL*«2 
RETURN 


WFP01S60 
WFP01S70 
WFPOlbHO 
WFP01590 
WFP01600 
taFPO  1 610 
WFP01620 
WFP01630 
WFPlI  1640 
WFPOlbbO 
wFPO 1660 
WFP01670 
WFP01680 
WFPO 1690 
WFPO 1/00 
WFP01710 
WFPO 1 720 
WFP01730 
WFP01740 
WFPO) 750 
WFPO 1760 
WFP01770 
WFPO 1780 
WFP01790 
WFP01800 
WFPOlblO 
WFP01820 
WFP01830 
WFP01840 
WFP01850 
WFP01860 
WFP01870 
WFPOlbbO 
WFP01890 
WFP01900 
WFP0191 0 
WFP01920 
WFPO 1930 
WFP01940 
WFP01950 
WFPC1960 
WFP01970 
WFP01980 
WFP01990 
WFP02000 
WFP020 1 0 
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SPECIAL  TEST  BLOCK  FOR  XL.GT, 
TEST  FOR  LEFT  INTERSECTION 

A=XD-XL 

IF  (A.lE.(-Ah))  U-SO 

(A.GT. (-Ah)  .ANO.A.lE.  <-AL)  )  14=80 
(A.GT.(-AL) .ANO.A.lE. (AL) )  14=70 
(A.GT .Al. ANO.A.lE. AH)  14=80 


AH 


IF 

IF 

IF 

If 


(A.GT.  AH)  14=9 


TESTS  INVOLVING  XMIn 


IF  (I4.EO.  700)  XMlN=(XD**2»XL**2-AL»*2)/(2.0*XD) 
IF  ( XMlN.LE.O.O)  13=1 
IF  (XMIN, GT .0.0)  13=2 


NUM=I4*I3 


determine  COMPONENT  SUM 


calculate  routing  argument 


I TST= i 

00  19  J=1 ,  8 

IF  (NUM.GF.MUMttER(J))  I T ST= I TS  ♦  1 
CONTINUE 

BRANCH  TO  >PPRQPRIATE  CODE  BLOCK 


GO  TO  (2.20, IS. 21 ,22.8.3) ,  ITST 

CODE  BLOCK  SECTION— 

XAREA=Pl*  (AH**<r-AL**2) 

RETURN 

AHANG= ANGLE  <  Ah.  XL), XL  )  *2.0 
XL ANG=AN(.Lt  (XL.XO.AH)  *2.0 
ASEGMT=SFGMT2(XL .XLANo) 

BSEGM1 =  Sf GMT2 ( AH , AhANu ) 

XLANG=ANGLE (XL.XO,AL)*2.0 
0PANG=AN0LE(XL,AL,XD) 

ALAN3=(XLAnG/2.0»OPANG) *2.0 
CSEGMT=Sf GMT2 ( xL.XLANG) 

DSEGMT  =Sf  GMT2 ( AL , Al ANG) 

XAREA=ASEGMT»8SEGMT»DSEGMT-CSFGMT-PI*AL**2 

RETURN 

AHANG=ANGLE (AH,XO.XL)*2.0 
XL ANG= ANGLE (XL ,X0, AH) *2.0 
ASEGMT  =  Sf GMT  2 ( XL, XL ANG) 
t)SEGMT=Sf  GMT 2  (  ah,  AHANG) 

XLANG  =  ANGLE(XL,X(),  AL)*2.0 
ALANG=AN  .,LE(AL,XD,XL)  *2.0 

(  segmt=si:gmt?( al .alang) 

i  :,EGMT=SEGMT2(XL.XLANG) 

XARE A=ASf  GMT»bSEGMT-CsEGMT-DSEGMT 
RF  TURN 

END 


MFP02020 
«IFP020  30 
MFP02040 
MFP020SO 
MFP020b0 
WFP02070 
WFP02080 
MFP02090 
MFP02100 
WFP021 10 
*F P02120 
WFP02130 
MFP02140 
*(FP021S0 
MFP02160 
WFP02170 
WFP021 80 
WFP02190 
RFP02200 
MFP02210 
WFP02220 
MFP02230 
MFP02240 
WFP022S0 
MFP02260 
MFP02270 
RFP02280 
MFP02290 
WFP02300 
WFP02310 
WFP02320 
WFP02330 
WFP02340 
WFP023S0 
MFP02360 
MFP02370 
WFP02380 
WFP02390 
MFP02400 
RFP02410 
RFP02420 
RFP02430 
RFP02440 
WF P024S0 
WFP02460 
KFP02470 
WFP02480 
MFP02490 
RFP02500 
WFP02S10 
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APPENDIX  IV 
SAMPLE  OUTPUT 
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^.00000 


AL=  3.00000 


AH=  10.00000 


XL  = 


XO  INTERSECTING  AREA 


0.00000 
.2<*000 
.  4  «  0  0  0 
.  72 000 
.*6000 
1.20000 

1  .44000 
1  .68000 
1.92000 
2.16000 
d.  40000 

2  .  64i)00 
2.88000 
3 • 1 200U 

3.36000 

3.60000 
3. *4000 
4.08000 
4.12000 
‘t  .56000 
*.*0000 
5.04000 
5  •  2  *  0  0  0 
6.52000 
5.76000 
6.00000 
o. 24000 
8.48000 
6. 72000 
6.96000 
7.20000 
7.44000 
7.68000 
7 .9? 000 
6.16000 
6.40000 
8.64000 
8.88000 
v. 12000 

4. 36000 

9.60000 
4.84000 

10.08000 
10.  J2uOO 
10.56000 
10.80000 
1 1 .04000 
1 1 .88000 
1  1  .6/90  00 
1 1 . 76000 
16.00000 


0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
. 37634 
1.12728 
1  .99964 
2.92944 
3.88356 
4.84280 
5.79372 
6.72599 
7.63074 
8.49974 
9.32468 
10.09655 
10.80490 
11.43653 
11.97297 
12.38292 
12.56637 
12.66637 
12.66637 
12.56637 
12.56637 
12.56637 
12.56637 
12.56637 
12.56637 
12.56637 
12.56637 
12.56637 
12.56637 
12.37909 
11. 84675 
11.15128 
10.3507? 
9.47901 
8.56029 
7.61 358 
6.65491 
5.69857 
4.75794 
3.84615 
2.97660 
2.16476 
1 .4276? 
. 78635 
.28254 
.00000 
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PK 

0.0000000 
0.0000000 
0.0000000 
0.0000000 
0.0000000 
.0013234 
.0039431 
.0069953 
.0102469 
.0135844 
.0169397 
.0202659 
.0235269 
.0266917 
.0297313 
.0326169 
.0353166 
.0377946 
.0400040 
.0418804 
.0433144 
.0439560 
.0439560 
.0439560 
.0439560 
.0439560 
.0439560 
.0439560 
.0439560 
.0439560 
.0439560 
. 0439560 
.0439560 
.0439560 
.0433010 
.0414389 
.0390062 
.0362059 
.0331567 
.0299431 
.026631b 
.0232783 
.0199331 
•  0 1 b6426 
.0134535 
.010412b 
.00  75721 
.0049937 
.0027576 
.0009883 
.0000000 


XL=  '3.00000  AL=  9.00000  AH=  10.00000 

«l)  intersecting  area  pk 


0.00000 

0.00000 

0.0000000 

.30000 

0.00000 

0.0000000 

.60000 

0.00000 

0.0000000 

.90000 

0.00000 

0.0000000 

1.20000 

0.00000 

0.0000000 

1.60000 

0.00000 

0.0000000 

1  .“*00  00 

0.00000 

0.0000000 

2. 1000U 

0.00000 

0.0000000 

2*40000 

0.0000° 

0.0000000 

2.  70000 

0.00000 

0.0000000 

J. 00000 

0.00000 

0.0000000 

3.30000 

0.00000 

0.0000000 

3.60000 

0.00000 

0.0000000 

3.9U000 

0.00000 

0.0000000 

4.20000 

.55168 

.0092457 

4.50000 

2. 10544 

.0352727 

4.80000 

4.11938 

.0690126 

5.10000 

6.24355 

.1045991 

5.40000 

7 , 50 1 bO 

.1256755 

5.70000 

8.36584 

.1401542 

6-00000 

9.00187 

.1508097 

6.30000 

9.47963 

.15881 36 

6.50000 

9.83843 

.1648247 

6.90000 

10.10337 

.1692633 

7.20000 

10.29168 

.1724181 

7.50000 

10.41576 

.1744967 

7.80000 

10.48483 

.1756540 

H. 10000 

10.50593 

.1760075 

H. 40000 

10.48450 

.1756483 

d. /0000 

10.42474 

.1746473 

9.00000 

10.32996 

.1730594 

9.30000 

10.20269 

.1709273 

4.60000 

10.04485 

.1682829 

9.90000 

9.85779 

.1651491 

10.20000 

9.64238 

.1615403 

10.60000 

9. 39901 

.1574631 

10.80000 

9.12761 

.1529162 

1 1 .10000 

8.82755 

.1478893 

1 1.40000 

8.49765 

.1423624 

11.70000 

8.13599 

.1363035 

12.00000 

7.73970 

.1296644 

12.30000 

7.30463 

.1223755 

12.60000 

6.82466 

. 1 143346 

12.90000 

6.29061 

.1053875 

13.20000 

5.68775 

.0952877 

1 3.60000 

4.98997 

.0835977 

1 3.H0000 

4.13992 

.0693567 

14.10000 

2.89877 

.0485635 

14.40000 

1.58538 

.0265601 

14. 70000 

.56312 

.0094340 

16.00000 

.00000 

.0000000 
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AH=  10.00000 


XL=  S. 00000  AL=  1.00000 


XU 

INTERSECTING  AREA 

PK 

0.00000 

75.39822 

.2424242 

.30000 

75.39822 

.2424242 

.60000 

75.39822 

.2424242 

.90000 

75.3982? 

.2424242 

1.20000 

75.39822 

.2424242 

1.50000 

75.39822 

.2424242 

i  .80000 

75.39822 

.2424242 

P.10000 

75.39822 

.2424242 

2.40000 

75.39822 

.2424242 

2.70000 

75.39822 

.2424242 

3.00000 

75.39822 

.2424242 

3.30000 

75.39822 

.2424242 

3.60000 

75.39822 

.2424242 

3.90000 

75.39822 

.2424242 

4.20000 

75.57813 

.2430027 

4.50000 

76.05929 

.2445497 

4.80000 

76.63642 

.2464054 

5.10000 

77.04684 

.2477250 

5.40000 

76.34  759 

.2454767 

8. 70000 

75.01119 

.2411798 

6.00000 

73.13757 

.2351557 

6. 30000 

70.75032 

.2274801 

6.60000 

68.18559 

.2192338 

6.90000 

65.48726 

.2105580 

7.20000 

62.68811 

.2015580 

7 .60000 

59.81375 

.1923162 

7.80000 

56.88488 

.1828992 

8.10000 

53.91870 

.1733622 

8. 40000 

50.92985 

.1637523 

8. 70000 

47.93106 

. 1541104 

9.00000 

44.93360 

.1444728 

9. 30000 

41 . 94765 

.1348722 

9.60000 

38.98257 

.1253388 

9.90000 

36.04712 

.1159005 

10.20000 

33.14960 

.1065843 

10.50000 

30.29808 

.0974159 

10.80000 

27.50045 

.0884208 

1 1 .10000 

24,76462 

.0796245 

11.40000 

22.09862 

.0710526 

11.70000 

19.51077 

.0627320 

12. 000 ''0 

17.00980 

.0546908 

12.30000 

14.60510 

.0469591 

12.60000 

12.30695 

.0395699 

12.90'  00 

10.12686 

.0325604 

13.20000 

8.07809 

.0259731 

13. ->0000 

6.17642 

.0198587 

1 3.80000 

4.44141 

.0142802 

14.10000 

2.89877 

.0093203 

14.4 0  0  0  0 

1.58538 

.0050974 

14. 70000 

.56312 

.0018106 

16.00000 

.00000 

.0000000 
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El 


■A  i 


.■■.niaOiir 


XL=  d. 00000 


AL=  7.00000 


10.00000 


X0 

INTERSECTING  AREA 

PK 

0.00000 

47.12389 

.2941176 

.36000 

47.123B9 

.2941176 

.72000 

47.12389 

.2941176 

1 .08000 

47.43221 

.2960420 

1.44000 

50.60648 

.3158538 

1.80000 

54.61001 

.3420895 

2.  16000 

58.68561 

.3662787 

2.52000 

60.22305 

.3758743 

2.88000 

60.86795 

.3798994 

3.24000 

61.03798 

.3809607 

1.60000 

60.90900 

.3801557 

3.06000 

60.57476 

.3780695 

4.32000 

60.09119 

.3750514 

4.68000 

59.49415 

.3713250 

5.04000 

58.80770 

.3670406 

5.40000 

58.04861 

.3623029 

5.76000 

57.22883 

.3571863 

6.12000 

56.35703 

.3517451 

6.48000 

55.43959 

.3460190 

6.84000 

54.48115 

.3400371 

7.20000 

53.48512 

.3338205 

7.56000 

52.45390 

.3273842 

7.92000 

51.38911 

.3207385 

8.26000 

50.29170 

.3138891 

8.64000 

49.16206 

.3068386 

9.00000 

48.00011 

.2995864 

9.36000 

46.80524 

.2921289 

9.72000 

45.57643 

.2844594 

10.08000 

44.31215 

.2765685 

10.44000 

43.01037 

.2684436 

10.80000 

41 .66648 

.2600684 

1 1 .16000 

40.28320 

.2514224 

11.62000 

38.85046 

.2424801 

1 1 .88000 

37.36514 

.2332097 

12.24000 

35.82086 

.2235712 

12.60000 

34.20945 

.2135138 

12.96000 

32.52032 

.2029714 

13.32000 

30.73930 

.1918554 

13.68000 

28.84664 

.1800425 

14.04000 

26.81302 

.1673500 

14.40000 

24.59041 

.1534779 

14.76000 

22.08238 

.1378243 

15.12000 

18.92856 

. 1181402 

15.48000 

15.54504 

.0970224 

15.84000 

12.37714 

.0772503 

16.20000 

9.44677 

.0589608 

16.56000 

6.781/3 

.0423273 

16.92000 

4.41918 

.0275817 

1 7.28000 

2.41325 

.0150620 

17.64000 

.85594 

.0053422 

13.00000 

.00000 

.0000000 
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1.00000 


AH= 


9.00000 


XL=  10.00000  AL= 


XO  INTERSECTING  AREA 


0.00000 
.38000 
.76000 
1.14000 
1.52000 
1.90000 
2.28000 
2.66000 
3.04000 
3.42000 
3.80000 
4.18000 
4.56000 
4.94000 
8.32000 
5.70000 
6.08000 
6.46000 
6.84000 
7.22000 
7.60000 
/. 98000 
8.36000 
8.74000 
9.12000 
9.50000 
9.88000 
10.26000 
10.64000 
11.02000 
1 1  .40000 
1  1.78000 
12.16000 
12.54000 
12.92000 
13.30000 
13.6H000 
14.06000 
14.44000 
14.82000 
15.200  'ft 
15.58000 
15.96000 
16.34  J 00 
16.72000 
17.10000 
17.48000 
l 7.86000 
18.24000 
18.62000 
19.00000 


251.32741 

251.32741 

251.32741 
250 . 44449 
245.79353 
239.98318 
233.69050 
227. 14894 
220.46681 
213.70291 
206.89277 
200.05973 
193.22024 
186.38658 
179.56842 
172.77372 
166.00931 
159.28123 
152.59502 
145.95583 
139.36856 
132.83798 
126.36876 
119.96554 
113.71375 
108.01243 
102.56294 
97.22090 
91.87607 
86.33947 
80.53349 
74.83394 
69.24679 
63.77836 
58.43537 
53.22502 
48.15502 
43.23375 
38.47036 
33.87492 
29.45865 
25.23417 
21.21592 
17. 42067 
13.86837 
10.58337 
7.59660 
4.94950 
2.70250 
.95841 
.00000 
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PK 

1.0000000 

1.0000000 

1.0000000 

.9964870 

.9779814 

.9548627 

.9298250 

.90379*9 

.8772095 

.8502969 

.8232002 

.7960124 

.7687989 

.7416086 

.7144800 

.6874448 

.6605300 

.6337599 

.6071563 

.5807398 

.5545299 

.5285455 

.5028053 

.4773277 

.4524526 

.4297678 

.4080850 

.3868297 

.3655633 

.3435338 

.3204326 

.2977548 

.2755242 

.2537660 

.2325070 

.2117756 

.1916027 

.1720216 

.1530687 

.1347840 

.1172123 

.1004036 

.0844155 

.0693147 

.0551805 

.0421099 

.0302259 

.0196934 

.0107529 

.0036134 

.0000000 


02/25/74  SCOPE  3.2.0  -  SCM  VER  130  43  08FE674ECPA=  1730 
08. 16.22. WFPNE73  AO 
OB. 16.22. WFPNEW.P6»T40.Cm60000. 

OB. 16.22. TASK (LOP  I  NO  * 88090904-9E M . SAH . 2295 ) 

0B.l6.22.FTN(B=E TNXAh.O) 

08.16.53.  3.485  CP  SECONDS  COMPILATION  TIME 

0  8 . 1 6. 54 .MAP (OFF ) 

08. 1 6. 54. RE W I  NO (OUTPUT) 

06. 16. 54. PRESET. 

08.16. 55  .FTNiXAR(LC=J  77777) 

08.17.00.  15200  CM 

08. 17. 02. STOP 

08. 17. 02. CP  004.255  SEC. 

08. 17. 02. PP  034.971  SEC. 
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